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Abstract. We report on the observation of hot water vapor (steam) in the inner AU of a young massive star located in the star¬ 
forming region IRAS 08576-4334. The water lines are detected in a medium resolution (R ~ 10,000) A"-band spectrum taken by 
the infrared spectrometer ISAAC mounted on the VLT-ANTU. The water vapor is at a mean temperature of 1565±510 K, cooler 
than the hot CO found in the same object, which is at = 1660 K and the column density is 1V(H20)= (2.5 ± 0.4) x 10 18 cm -2 . 
The profile of both TEO and CO lines is best reproduced by the emission from a Keplerian disk. To interpret the data, we also 
investigate the formation of molecules and especially CO and water vapor in the inner hot and dense part of disks around young 
high mass stars using a pseudo time-dependent gas-phase chemical model. Molecules are rapidly photodissociated but this 
destruction is compensated by an efficient formation due to fast neutral-neutral reactions. The ability of CO molecules to self¬ 
shield significantly enhances its abundance. Water molecules are sufficiently abundant to be detectable. The observed H 2 0/CO 
ratio is reproduced by gas at 1600 K and an enhanced UV field over gas density ratio /uv/rcH=l(L 4 -10~ 6 . The simulations 
support the presence of CO and H 2 0 molecules in the inner disks around young massive stars despite the strong UV radiation 
and show that the OH radical plays an essential role in hot gas chemistry. 
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1. Introduction 

Disks are found around many low-mass young stars and are 
natural byproducts of the star formation process (e.g.. Waters 
& Waelkens 1998; Natta et al. 2000 for the Herbig Ae stars 
and Mundy et al. 2000 for the T Tauri stars). Unfortunately, 
observations of disks around their higher counterparts remain 
inconclusive (Churchwell 2002). The presence or absence of a 
gaseous disk may help understanding the detailed physics of 
high-mass star formation. One of the main difficulties in the 
formation of massive stars by direct molecular core collapse 
is that the radiation pressure of a luminous massive protostar 
will halt the accretion of matter, ultimately limiting the final 
mass of the star (Wolfire & Cassinelli 1987). Two major sce- 
narii have been proposed to overcome this shortcoming. In the 
first one, massive stars can form by coagulation of two lower 
mass objects in a dense star cluster (Bonnell & Bate 2002). 
In the second scenario, a massive circumstellar disk will allow 
a significant amount of matter to accrete to the star (Yorke & 
Sonnhalter 2002; Jijina & Adams 1996). Evidence for the pres- 
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ence of a dense gaseous disk close to the star NGC2024-IRS2 
were presented by Lenorzer et al. (2004). The main argument is 
that the infrared excess observed in NGC2024-IRS2 can be best 
explained by the emission from a dense gaseous disk extending 
up to 0.6 AU. The knowledge of the structure and composition 
of the inner disk is essential for our understanding of the accre¬ 
tion of matter onto the star. At the high densities and tempera¬ 
tures (100-2000 K) characteristic of disks at 0.1-5 AU around 
young stars, molecules are expected to be sufficiently excited 
to produce a rich ro-vibrational spectrum in the near- and mid- 
infrared. CO overtone emission, which is most likely emitted 
by a disk, has been detected and analyzed by Scoville et al. 
(1983) toward Orion BN object and by Chandler et al. (1995) 
toward NGC2024-IRS. Disk models have been successful in 
modeling CO bandhead emission in low-mass as well as high- 
mass young stellar objects (Najita et al. 1996; Chandler et al. 
1995; Kraus et al. 2000). Cooler CO fundamental emission, 
which probes the planet forming region of disks, has been de¬ 
tected in a large number of low-mass pre-main-sequence stars 
(Blake & Boogert 2004; Najita et al. 2003; Brittain et al. 2003). 

Recent observations of CO overtone emission from young 
early type stars further support the disk hypothesis (Bik & Thi 
2004; Blum et al. 2005). The presence of CO molecules raises 
the question of their formation and/or survival in the harsh 
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environment in the close vicinity of massive stars. Molecules 
other than CO are also expected to exist in the inner few AU 
of disks. Carbon, nitrogen and oxygen are the most abundant 
elements after Hydrogen and Helium. At high temperatures 
and low pressures the dominant equilibrium species of these 
elements are CO, N 2 , and H 2 O; at intermediate temperatures, 
CH 4 , N 2 and H 2 O; and at low temperatures and high pressures, 
CH 4 , NH 3 , and H 2 O (Lewis & Prinn 1980). In all cases wa¬ 
ter vapor is the second molecular reservoir of oxygen after CO 
in the gas phase at a temperature between its dissociation tem¬ 
perature (~2500 K) and its condensation temperature (~ 100- 
150 K). Basically, all gas-phase oxygen locked into H 2 O does 
not participate to further reactions (e.g., Kaufman & Neufeld 
1996; Elitzur 1979). Chemical models of protoplanetary disks 
abound in the literature for low-mass stars (e.g., Aikawa et al. 
1999; Markwick et al. 2002; Semenov et al. 2004), high-mass 
stars (e.g., Nguyen et al. 2002) or debris-disks (e.g., Kamp et al. 
2003) although most studies were limited to the region located 
beyond a few AU. Nevertheless, several studies focus on the 
thermal-chemical structure of the inner disks around young low 
mass stars (Glassgold et al. 2004). 

Fundamental ro-vibrational and pure-rotational water lines 
have been observed in the envelope of high-mass young stellar 
objects using the Short- and Long-Wavelength Spectrometer 
on board the Infrared Space Observatory (e.g., Boonman et al. 
2003; Wright et al. 2000; Helmich et al. 1996). They find that 
freeze-out of water onto grains predominates in the cold part 
of the envelope while the high abundance of water in the warm 
part can only be explained by sublimation of water molecules. 

Water emission lines have also been observed in the sub¬ 
millimeter domain with the Submillimeter Wave Astronomy 
Satellite (SWAS) but the beam was too large to allow the study 
of small objects like the disks around high mass young stellar 
objects. Surprisingly, SWAS data show that gaseous H 2 O and 
O 2 are not the dominant carriers of elemental oxygen in cold 
molecular clouds; the oxygen remains in atomic form, provid¬ 
ing strong constraints on chemical models (Melnick et al. 2000; 
Bergin et al. 2000). 

Highly excited water lines, which are absent in our atmo¬ 
sphere, are detected from the ground in the atmosphere of M 
dwarfs (see Allard et al. 2001 and references therein). Sunspots 
(Zobov et al. 2000) and in the disk around the low-mass young 
star S VS 13 (Carr et al. 2004). In the latter object, CO bandhead 
emission has also been detected. 

In this paper, we model the hot water emission 
lines detected in the disk around the massive young star 
08576nr292 (IRAS 08576-4334). The young stellar object 
(YSO) 08576nr292 is one of the infrared sources detected in a 
survey of Hii regions by Kaper et al. (2005). We also investigate 
the condition for the formation and/or survival of molecules 
(mainly CO, H 2 O) in a gaseous disk around a massive young 
star. The inner disk around a young massive star is probably 
too hot for dust grains to survive, such that no H 2 formation on 
the grain surface can occur. Moreover, far ultraviolet photons 
are no longer absorbed by dust so that photodissociation oc¬ 
curs at high column densities. In addition, the disk temperature 
and density are high, making the chemistry different compared 
to other environments. For example, endothermic reactions or 


reactions with an activation energy, mostly the neutral-neutral 
reactions, and three-body reactions may dominate the chemical 
network over the classical ion-molecule reactions. The chem¬ 
istry scheme is probably closer to that in shocks or protostellar 
winds than that of a cold quiescent molecular cloud. 

The paper is organized as follows: we first present the ob¬ 
servations and the fit to the K-band spectrum in Sect. 2. Then 
we detail the chemical models in Sect. 3. In Sect. 4, we present 
and discuss the results of the numerical runs. We finally con¬ 
clude in Sect. 5. 

2. Observations 

2.1. The high-mass young star 08576nr292 

IRAS 08576-4334, identified by Liseau et al. (1992) as IRS 
34 in their survey of IRAS sources in the Vela Molecular 
Cloud, is a high-mass star forming region with total bolomet- 
ric luminosity log (L/L Q )= 4.35. The distance to the cloud has 
been estimated to be around 700 pc. The extended Hu region 
IRAS 08576-4334 includes a radio continuum source known 
as G265.151+01.454 (Walsh et al. 1999). 08576nr292 is one 
of the near-IR sources detected in the J- and K-band image 
of IRAS 08576-4334 by Kaper et al. (2005) and discussed in 
Bik et al. (2004). The properties of 08576nr292 are listed in 
Table 1. After correction by 12 magnitude of extinction, as de¬ 
rived by Bik et al. (2004), the Jo - Kq color of 08576nr292 is 
small (0.4). Because the intrinsic color of hot stars is —0.2, the 
near-IR excess from 08576nr292 is (J - K)= 0.6, suggesting no 
strong hot dust emission. 


08576nr292 


Ra (72000) 

08:59:21.6 

Dec (72000) 

-43:45:31.6 

K (Mag) 

9.4 ± 0.01 

7 - K (Mag) 

4.7 ± 0.09 

Am 

12 

K q (Mag) 

-3.1 

J 0 - K 0 (Mag) 

0.4 

7(kpc) 

0.7 

M„(M 0 ) 

6 


Table 1. Properties of 08576nr292. The position and observed magni¬ 
tude and color is taken from Kaper et al. (2005). The extinction esti¬ 
mate and the derived unreddened magnitude and color are taken from 
Bik et al. (2004). 


2.2. Observation & data reduction 

The spectrum of 08576nr292, identified in the line of sight to 
IRAS 08576-4334, was observed with ISAAC at the Very Large 
Telescope ANTU in 2002 at resolving power of 10,000 (slit 
width of 0.3"). The near-infrared data were reduced in a stan¬ 
dard way using a combination of the dedicated software pack¬ 
age Eclipse and IRAF. We made use of Hatfield taken during 
the night and arc frames were taken by the European Southern 
Observatory (ESO) staff during the day. 
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Standard stars of spectral type A observed at similar air- 
mass were used to correct for telluric OH emission and ab¬ 
sorption. The wavelength calibration was achieved by using arc 
lines in combination with atmospheric OH lines. The details of 
the observations and datareduction are presented in Bik et al. 
(2004). The chosen setting covers the spectral range from 2.280 
till 2.400 micron so that the four first CO bandhead (Av=2, 
where v is the vibrational quantum number) can be observed. 
At resolving power 10,000 (Av - 30 km s -1 ), the profile of the 
lines is not resolved and a few lines may even be blended. 



Wavelength (micron) 


Fig. 1. The lower panel shows the best fit to the CO first bandhead to¬ 
ward 08576nr292. In the upper part the observed spectrum is shown. 
The spectrum is continuum subtracted and normalized to the peak 
of the 2-0 bandhead. The best fit of the CO-bandhead is plotted in 
the lower part. The parameters of the fit are discussed in Bik & Thi 
(2004). Features shortward of 2.2935 pm are shown in detail in the 
upper panel. 


2.3. Observed water lines 

The spectrum shows four prominent CO bandhead emissions 
which are modeled in Bik & Thi (2004). The fit to the first 
CO bandhead is shown in the lower panel of Fig. 1. The upper 



Fig. 2. Best fit to the hot water lines with the slab model. The up¬ 
per spectrum is a blowup of the spectrum showed in Fig. 1 between 
2.278 and 2.292 pm, normalized to the maximum strength at ^2.280 
pm. The middle one represents the best fit while the lower curve is 
the synthetic spectrum before degradation at the resolution of the ob¬ 
served spectrum. The characteristic of the strong detected lines, which 
are labeled 1 to 10, are given in Table 2. The Ho 3-2 S(l) emission line 
is located at position 6. 


panel displays in detail the spectrum shortward of A < 2.2935 
pm. The CO bandhead is best fitted by gas at 1660 K located 
in a Keplerian rotating disk extending from 0.2 till 3.6 AU and 
viewed with an angle of 20-30° with respect to the disk rotation 
axis. The column density of CO is relatively large (A(CO)^3.9 
x 10 21 cm -3 ). 

Several emission features can be seen blueward of the 
first CO bandhead (i.e. A < 2.2935 pm), which we attribute 
to water vapor at T>1000 K. The emission features were 
checked against atmospheric transmission spectra to ensure 
no artifact is present. Using spectroscopic parameters drawn 
from the compilations HITEMP (Rothman et al. 2003) and 
that of Tennyson et al. (2001), most lines (often blended) in 
the 2.280 - 2.2935 pm region can be assigned to hot water 
emission lines. HITEMP (High-Temperature molecular spec¬ 
troscopic database) contains water vapor spectroscopic param¬ 
eters of lines which have strength greater than 3 x 1 0 27 
cm '/(molecule cm -2 ) at 1000 K while the energy of the lev¬ 
els in Tennyson et al. (2001) were carefully assessed using 
energy levels derived from sophisticated variational methods. 
The HITEMP database contains in total 1816 water lines be¬ 
tween 2.278 and 2.292 pm which are sufficiently strong at 
r ga s=1500 K, but only 10 lines have strength greater than 1 0 2 5 
cm '/(molecule cm -2 ), which is likely the threshold for visual 
detection in the observed spectrum. The detected strong wa¬ 
ter lines and their spectroscopic characteristics are summarized 
in Table 2. The wavelength of three of the features (lines 8, 
9 and 10) matches those of hot water lines detected by Carr 
et al. (2004) in the inner disk around the low mass young star 
SVS 13. The spectrum of SVS 13 obtained by Carr et al. (2004) 
is at a higher resolution but covers a much smaller spectral 
range. 
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Table 2. Strongest detected water emission lines. Various spectroscopic characteristics of the detected lines are also provided. The last column 
lists estimates of the critical density ;i crit assuming r gas =1000 K and <x ^10 -15 cm -2 . 


Line 

Nuclear state 

Transitionfu—>1) 

(VH 2 V 3 ) J K a K c - (vjv'v') J' K' a K' c 

Wavelength 

(micron) 

^upper 

(cm -1 ) 

Slower 

(cm -1 ) 

A u i 

(s -1 ) 

^crit 

(cm -3 ) 

1 . 

ortho-Hi O 

(001) 15 6 10-(000) 14 4 11 

2.28020 

7131.63 

2746.02 

9.50 10 -2 

9.50 10 8 

2 . 

ortho-Hi O 

(001) 15 7 9-(000) 14 5 10 

2.28076°’ 

7302.72 

2918.24 

8.88 10 -2 

8.88 10 s 

3. 

ortho-HiO 

(001) 16 5 11 -(000) 15 3 12 

2.28093°’ 

7464.48 

3080.18 

8.27 10 -2 

8.27 10 8 

4. 

ortho-HiO 

(001) 14 8 6 -(000) 13 6 7 

2.28449° ’ 

7133.78 

2756.41 

6.56 10 -2 

6.56 10 8 

5. 

ortho-HiO 

(100) 18 7 12-(000) 17 4 13 

2.28462°' 

8394.96 

4017.90 

2.00 10 -1 

2.00 10 " 

6 . 

ortho-HiO 

(001) 12 9 3-(000) 11 7 4 

2.28694 

6694.58 

2321.90 

1.88 10 -2 

1.88 10 8 

7. 

ortho-HiO 

(001) 14 5 10-(000) 13 3 11 

2.28744 

6619.79 

2248.06 

6.24 10 -2 

6.24 10 8 

8 . 

ortho-HiO 

(001) 13 8 6 -(000) 12 6 7 

2.29036 

6799.95 

2433.80 

4.65 10 -2 

4.65 10 8 

9. 

ortho-HiO 

(001) 13 4 10-(000) 12 2 11 

2.29134 

6139.03 

1774.75 

3.20 10 -2 

3.20 10 8 

10 . 

ortho-HiO 

(001) 18 7 11 -(000) 17 5 12 

2.29180 

8537.65 

4174.38 

3.21 10 -1 

3.21 10" 


Note 1 Blended lines. 


It should be noticed that all 1816 lines are included in the 
synthetic spectra described in Sect. 2.4. The position of the 
lines is indicated in Fig. 2 and Fig. 4. 

All strong detected lines arise from the first vibrational ex¬ 
cited symmetric ( 100 ) and antisymmetric ( 001 ) stretch level 
of ortho-H 20 . It is known that there is a near coincidence of 
the frequencies of those two modes. The rotational levels are 
high (12 < J < 18) as expected since only high level lines 
are not present in our atmosphere and can be readily detected 
with ground telescopes. At the resolving power of 10,000 (or 
resolution Av - 30 km s 1 ), the lines are unresolved and few 
of them blended. The observed lines arise from levels with en¬ 
ergy from 5000-9000 K. The detected (001)-(000) transitions 
have A K a = -2 and A K c = +1 and are therefore of o-type. The 
(lOO)-(OOO) transition is of />type (A K lx - -3 and A K c = +1) 
(see Table 4). 

Water lines are clearly seen in this object because they are 
not blended. The profile of emission lines coming from the 
inner part (< a few AU) of a disk in Keplerian rotation can 
have FWHM up to =* 200-250 km s ' 1 (or AT = 0.0015 /urn) 
if viewed edge-on. In comparison water lines wavelength po¬ 
sitions are often separated by less than 65 km s ~ 1 (or 0.0005 
yum). Therefore the rotational broadening of all lines produces 
a blend of features that can be mistaken for a continuum. 


2.4. Spectral modeling 

In this section we apply two different models to fit the observed 
spectrum. The first model is a simplified version of the disk 
model applied to the CO bandheads of this object by Bik & Thi 
(2004). In this model we assume that the gas is originating from 
an isothermal slab of gas. The second model which is used to 
fit the data is similar to the model used to fit the CO bandheads. 
The gas is located in a Keplerian rotating, isothermal disk. First 
we describe the ro-vibrational transitions of water vapor and 
the transition selection rules. 


2.4.1. H 2 0 ro-vibrational transitions 

The structure of the water molecule is described in molecular 
physics textbooks (e.g., Bernath 1995). As a triatomic asym¬ 
metric (near-oblate) top molecule, water has three vibrational 
modes and three rotational modes. The vibrational levels are 
characterized by three vibrational quantum numbers Vi (sym¬ 
metric stretching mode vi), v 2 (bending mode v 2 ) and V 3 (asym¬ 
metric stretching mode V3). The energy of the lowest vibra¬ 
tional levels is given in Table 3. The rotational levels are labeled 
using the standard asymmetric top notation Jk„k c , where J is 
the main rotational quantum number, K a is the prolate quan¬ 
tum number (projection of J along the a axis) and K c the oblate 
quantum number (projection of J along the c axis). Notice 
that K a + K c = J or / + 1. The rotational statistical weight 
is (gj = 2 J + 1). The total degeneracy is g - gjgj. 

Like H 2 , H 2 0 has two hydrogen atoms each having a nu¬ 
clear spin 1/2, which combine to a singlet (para-H 2 0) and 
triplet (ortho-H 2 0). Transitions between para and ortho states 
are highly forbidden. Ortho states have a nuclear statistical 
weight gi of 3 while para states have gj equal to 1. The sum 
K a + K, + i ’3 is odd for ortho states and even for para states. 

The selection rules for asymmetric tops depend on the com¬ 
ponents (yu a , /.(b, fj c ) of the permanent dipole moment vector, 
along the a, b and c principal molecular axes. The vibrational 
selection rules are Av = 1 for the strong fundamental transi¬ 
tions and Av = 2, 3, 4, ... for the weaker overtone transitions. 
The rotational selection rules can be divided into three general 
cases given in Table 4 in addition to AJ — 0, ±1. 

2.4.2. Slab model 

For the spectral modeling of the water lines we employ a 
method similar to the fit of the CO bandhead emission. The 
CO bandhead emission is best modeled by a Keplerian disk 
seen at low inclination with respect to the rotation axis. First 
we further simplify the model and assume that the water va¬ 
por lines arise from a isothermal slab of gas at temperature r gas 
with a column density A'(H 2 0). The remaining free parameter 
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Table 3. Lower energy levels of vibration of water molecule in cnr 1 
with ground energy level at 0 cnT 1 . The zero-point energy is 4500 
cm -1 . To convert the energy scale into K, multiply by 1.438768. 


Symmetry 

Vibrational level 
(viv 2 v 3 ) 

Energy 

(cm -1 ) 

Level type 

Ai 

000 

0000.00 

ground 

Ai 

010 

1594.75 

fundamental 

Ai 

020 

3151.63 

overtone 

Ai 

100 

3657.05 

fundamental 

b 2 

001 

3755.92 

fundamental 

Ai 

030 

4666.79 

overtone 

Ai 

110 

5234.97 

combination 

b 2 

011 

5331.27 

combination 

Ai 

040 

6134.01 

overtone 

Ai 

120 

6775.09 

combination 

b 2 

021 

6871.52 

combination 

Ai 

200 

7201.54 

overtone 

b 2 

101 

7249.81 

combination 

Ai 

002 

7445.04 

overtone 


Table 4. Selection rules for the ro-vibration transitions in addition to 
A J = 0, ±1. The transitions in brackets are weaker than the main 
transitions. For ro-vibrational transitions, the parity of K a + K c + V 3 has 
to be conserved. Transitions of type a only occur for oblate molecules 
whereas transitions of type c are only found for prolate molecules. 
Both oblate and prolate molecules can have b-type transitions. 


Transition 

Criterion 

\K, 

AK C 

a-type 

3* 0 

0, ±2 (±4, ...) 

±1 (±3, ±5, ...) 

fe-type 

Hb* 0 

±1,±3 (±5, ...) 

±1,±3 (±5, ...) 

c-type 

* 0 

±1 (±3, ±5, ...) 

0, ±2 (±4, ...) 


is the turbulent velocity width dv. Values for dv range from 5 
to 25 km s 1 . One should keep in mind that the profiles are 
unresolved and thus dv is an ill-constrained parameter. The tur¬ 
bulent profile is a Gaussian function. For dv >7 km s 1 , the 
gas is highly turbulent. The transition data are taken from the 
compilation HITEMP. The level population is assumed to be at 
local thermodynamical equilibrium (LTE). 

An order of magnitude check of this last assumption can 
be performed. The collisional de-excitation coefficients for 
the transition u—>1 is given by y u 1 ^ cr u iv re i, where cr ul is 
the de-excitation cross section and v re i = J%kT gas lmn\i 2 o = 
10 5 x/Tg as /1000 cm s 1 is the mean relative velocity of FLO 
molecules. In the absence of published cross-sections for the 
transitions considered here, we assume that <x u i is close to the 
geometrical cross section cr ~ 10 1 s cm 2 of the molecule. 
The rates of spontaneous emission can be estimated from the 
Einstein coefficients A u i, which range from 10 5 s 1 till 3 x 
10 2 s _1 in the V-haud region. The critical densities are approx¬ 
imately n crit = A u i/y u i ^ (10 5 -3xl0 12 )x(7’ gas /1000)~ 1/2 cm -3 . 
The number density in the inner part of disks is greater than 
10 10 cm -3 from the fit to the CO bandhead by a LTE model. 
Therefore only the strongest lines in the K -band are in non- 
LTE. More specifically, the lines detected here have moder¬ 


ately strong Einstein coefficients (10“ 2 —10 1 s -1 ) and by conse¬ 
quence relatively low critical densities n cnl which are reported 
in the last column of Table 2. 



Fig. 3. Typical line optical depth for the hot water lines with dv= 5 km 
s -1 , 7^=1500 K and A(CO)=3 x 10 18 cnr 2 . 


The best model fit to the data was obtained by minimizing 
the x 1 using a Simplex (Press et al. 1992) and a genetic algo¬ 
rithm methods ( Pikaia , Charbonneau 1995). The results of the 
two methods were cross-checked to ensure that no local min¬ 
ima were reached. The best fit using a slab model gives FLO 
gas at 1550 + 510 K and column density of (2.9 ± 0.4) x 10 18 
cnr 2 . The best fit parameters are gathered in Table 5. Although 
the errors on the gas temperature and column density are dom¬ 
inated by systematics, we estimated the uncertainties in the 
two parameters following the method described in Bevington 
& Robinson (2003). The derived errors represent here an in¬ 
crease of x 2 by 1 from its value at the minimum. Noteworthy, 
the statistical significance of the errors is difficult to asses. The 
column density of water is 3 orders of magnitude lower than 
CO while the temperature is only slightly lower for water than 
CO. The intrinsic turbulent width dv is an ill-constrained pa¬ 
rameter although a value of 5 km s 1 has been found to provide 
the best results. On the other hand the column density and gas 
temperature are well constrained. 

2.4.3. Disk model 

The second model includes broadening effects caused by a disk 
viewed with an angle of 27°, as found by fitting the CO band- 
head. The slab model spectrum is convolved with a profile 
which takes into account the Doppler shifts caused by hot wa¬ 
ter located in a Keplerian disk. The disk extends from R mm till 
R max which are two parameters of the model. The inclination 
is assumed equal to the value obtained by fitting the first CO 
bandhead (7=27 °). Line emission from a Keplerian disk has a 
double-peak profile. The simple convolution procedure is jus¬ 
tified for FLO column densities in the range 10 18 -10 2 ° cm -2 
where the optical depth is t <<1 (see Fig. 3). The central star 
is assumed to have mass of 6 M Q (Bik et al. 2004). 
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2.280 2.282 2.284 2.286 2.288 2.290 2.292 

Wavelength (micron) 


Fig. 4. Best fit to the water lines toward 08576nr292 using the disk model. The upper curve is the observed normalized and continuum subtracted 
spectrum. The middle one represents the best fit while the lower curve is the synthetic spectrum before degradation at the resolution of the 
observed spectrum. Notice the double-peak profile for lines which are not blended. 

Table 5. Parameters derived from the fit by the slab and disk model. The parameters for the CO bandhead fitting are derived in Bik & Thi 
(2004). The accuracy in the temperature and column density estimates is limited by systematic instead of statistical errors (see text for a fuller 
discussion). 


Model 

Species 

i 

n 

dv 

(km s~‘) 

7gas 

(K) 

N 

(crrr 2 ) 

^min 

(AU) 

*MTiaX 

(AU) 

Disk 

CO 

27 

5 

1662 

3.9 x 10 21 

0.2 

3.6 

Slab 

h 2 o 

n. a. 

25 

1580 ±510 

(2.9 ± 0.4) x 10 18 

n. a. 

n. a. 

Disk 

h 2 o 

27 

5 

1550 ±510 

(2.1 ± 0.4) x 10 18 

2.0 ± 1 

4.0 ± 1 


The dust sublimation radius, i.e. the radius where dust 
grains attain 1500 K, should provide a upper boundary to outer 
radius 7? max because beyond this radius dust opacity will dom¬ 
inate over the gas opacity and water nor CO lines can be ob¬ 
served. A total column density of gas of ~10 25 cirT 2 , derived 
by multiplying the CO column density by 10 4 results in a dust 
opacity of t ~ 500 in the AT-band using the extinction law 
of Weingartner & Draine (2001), which assumes a standard 
gas-to-dust mass ratio, and Rv-3.1. A CO abundance of 10 4 is 


reached for a large range of physical conditions (see Section 3). 
However, gas located in the upper atmosphere of the disk can 
be warmer than the mid-plane. The total luminosity toward 
08576nr292 is estimated to be around 10 3 L 0 for a star of mass 
- 6.0 M a . A crude estimate to the dust sublimation radius is 
r& = 0.0346 y[LJL ( ~ AU ^ 1 AU, where a dust sublimation 
temperature of 1500 K is assumed. 

The results of the fit by a disk model is displayed in Fig. 4. 
The quality of the data and the resolution do not permit an 












































W.-F. Thi and A. Bik: Steam in the circumstellar disk around a massive YSO 


7 


accurate constraint on the parameters such as R nnn and /? max . 
The gas temperature and column density are however better 
bracketed because these two parameters are set by the relative 
strength between lines at LTE. The best fit parameters are pro¬ 
vided in Table 5 together with that derived from fitting the first 
CO bandhead. 

2.4.4. Discussion on the spectrum fitting 

The slab and disk models give similar values for the water gas 
temperature and column density. This closeness can be ascribed 
to the fact that the lines are optically thin. The intensity of opti¬ 
cally thin emission lines is not sensitive to the actual geometry 
of the object. A mean value between the two models for the gas 
and column density is adopted in this discussion. The resulting 
error on the gas temperature is the scatter between the two tem¬ 
peratures and is much smaller than in the formal error analysis 
(15 K instead of 510 K). The water has a mean temperature of 
1565 + 15 K and a column density of (2.5 + 0.4) x 10 18 cm -2 . 
Both models give relatively similar values for the temperature 
and column density. This proximity can be ascribed to the low 
optical thickness of the water lines. 

The fits using either a slab or a disk model match relatively 
well the moderate signal-to-noise ratio spectrum except that 
both fail to reproduce the strength of the peak at 2.287 /mi 
(line 6) relative to that of at 2.2875 (line 7). The upper level 
energy of both lines are similar such that non-LTE effect can 
be excluded (see Table 2). Moreover the Einstein A coefficient 
is slightly lower for line 6 than 7. Therefore line 7 should ap¬ 
pear stronger in most situations. The most likely explanation to 
this discrepancy is that the observed emission at 2.287 //m is 
a combination of hot water and FE 3-2 S(l) emission. Indeed, 
a strong FE 1-0 S(l) line has been detected in 08576nr292 at 
2.218 /tm (Bik et al. 2004). The origin of the FE lines is not 
well determined and a large contribution from the extended Hu 
region is probable. Alternatively, errors and uncertainties in the 
Einstein-A coefficients may exist and can be sufficiently large 
to explain the bad fit around 2.287 /urn. 

The fit by a slab requires a higher turbulent width (dv =25 
km s 1 ) than the disk model (civ=5 km s 1 ). This difference 
stems from the difficulty to constrain dv with unresolved lines. 
The unlikely high value of dv for the slab model and the si¬ 
multaneous goodness of the fit to the CO and water lines with 
the Keplerian disk model favor the latter interpretation for the 
geometry of the hot CO and water emitting region. Synthetic 
spectra obtained with a Keplerian disk and supersonic values 
for dv give unsatisfactory fits. Therefore we can rule out a tur¬ 
bulent Keplerian disk. 

The water vapor column density of (2.5 + 0.4) x 10 18 
cm~ 2 and the gas temperature 7 v;lpor = 1565 K are close to 
the values found around the low-mass young star SVS 13 by 
Carr et al. (2004) which is a priori unexpected since SVS 13 
and 08576nr292 differ in mass. This similarity stems from that 
fact that the FEO emission comes from a region very close to 
SVS 13. Carr et al. (2004) find that the emitting area extents out 
to 0.3 AU while the water vapor is emitted from 2 to 4 AU in 


the disk around 08576nr292. The hot gas is located at a smaller 
distance from the low-mass star than from the massive star. 


The water vapor temperature is about 100-200 K lower 
than that of the hot CO, in accordance with the findings of Carr 
et al. (2004) for SVS 13, although the hot CO is much warmer 
around SVS 13. This lower temperature is also consistent with 
a large inner and outer radius (Table 5). Interestingly, the water 
temperature is close to the maximum dust sublimation tempera¬ 
ture of 1500 K and the outer radius R m i n is much larger than the 
dust sublimation radius of 1 AU, though the effect of large dust 
opacity in the observed spectrum is absent. In principle the low 
contrast between the line and the continuum should have pre¬ 
vented any possible detection of the weak optically thin water 
emission lines in the presence of dust. One possible explana¬ 
tion is the lack of dust grains close to the star at 1300-1500 K. 
The large luminosity of a B5 stars can push the dust grains well 
beyond the dust sublimation radius through radiation pressure. 
Scaling the figures computed by Saija et al. (2003) to the radius, 
mass and effective temperature of a B5 star, we obtain values 
for the ratios of radiation pressure force over gravity /3 which 
are much larger than 10 for all types of grain compositions. We 
have not taken into account the viscous drag generated by the 
friction between the dust grains and the ambient gas which may 
be important in the inner dense disk. The absence of hot dust is 
also compatible with the relatively small near-infrared excess 
E(J - K). As a result of the separation between the gas and the 
dust, the hot gas at T -1500 K would be deprived of elemental 
oxygen which comes from the sublimation of silicate grains. 
A lower oxygen abundance in the gas phase may significant 
affects the chemistry. We postpone the discussion on the chem¬ 
istry to Section 4.5. Another possible explanation for the low 
continuum optical depth resides in the difference in dust grain 
lifetime. The maximum grain lifetime at a certain temperature 
is proportional to its radius. Lamy (1974) estimated that the 
maximum lifetime of a 1 /mi radius silicate grain at 1300 K 
is -0.6 year and ten times shorter for a 0.1 /mi radius grain. 
Therefore, we can envision that only the larger grains in the 
dust size distribution can survive in the region where T =1300- 
1500 K, which results in a lower dust opacity. Alternatively, 
if the thermal coupling between the gas and the dust is ineffi¬ 
cient, especially in the upper atmosphere of disks, the gas can 
be hotter than the dust (Glassgold et al. 2004). In addition, the 
disk upper atmosphere can be relatively dust-free due to dust 
settling. This latter possibility may be at work because the disk 
model with a subsonic turbulent velocity is favored. Finally, it 
should be kept in mind that a combination of dust settling and 
size segregation is possible. 


We have demonstrated the presence of hot CO and FEO 
in the inner 4 AU of a Keplerian disk orbiting a B5 star. In 
the next section, a detailed chemical model is built in order to 
establish whether molecules can form and survive in the harsh 
environment around a hot star. 
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3. Chemical modeling 

3.1. Pure gaseous disk 

The presence of relatively large amount of HiO and CO in the 
inner disk around young early type stars raises the question of 
their origin. 

Early B stars are the main contributors to the mean inter¬ 
stellar radiation field. A critical part of the UV spectrum is the 
912-1000 A range used for the computation of FE and CO pho¬ 
todissociation rates. We modeled the chemistry that occurs in 
the HI regions using a pseudo-time dependent chemical code, 
i.e. all photons below 912 A have been absorbed. The main 
parameters of the models are the gas temperature r gas and den¬ 
sity «h, the strength of the UV field and the abundances of the 
species at time t- 0. The temperature (7' gas = 1500-4000 K) and 
density range f/? H = 10 1 f> — 10 14 cm -3 ) considered in our model¬ 
ing correspond to the values derived by Lenorzer et al. (2004) 
when fitting the spectral energy distribution of NGC2024/IRS2. 

The radiation field experienced by a parcel of gas is domi¬ 
nated by the emission from the central star diluted by the dis¬ 
tance to the star and by self-shielding. Because of the proximity 
to the central star, the disk can be considered flat and the radi¬ 
ation flux received by the disk surface at distance R from the 
star of radius R, with effective temperature T, is: 

fAtar = \(rTt (y) Ho = crTtW dlsk , (1) 

with 




( 2 ) 


Ho is the cosine of the angle between the incoming stellar radia¬ 
tion and the normal to the surface of the disk. For a typical B5 V 
star, R, — 3.9 R Q , therefore the dilution factor Wdisk at 2 AU is 
Wdisk ~ 7.5 x 1(T 7 compared to Wism=10~ 14 for the mean in¬ 
terstellar UV field. 

The photodissociation rates in the UMIST database are 
multiplied by a factor 7jjv = Wdi s k/wiSM=7.5 x 10 7 . The spec¬ 
trum of the mean interstellar UV field is close to that of a BOV 
star and it is assumed here that a unique enhancement factor 
suffices to simulate the UV photodissociation. We explore the 
effect of an enhanced UV field by factor /uv=10 6 -10 9 because 
of the uncertainty in the actual spectral type of 08576nr292 
(Bik et al. 2004) and in the inner and outer radius of the emit¬ 
ting area. In photodominated chemistry, an important param¬ 
eter is the ratio 7 uv/hh=10~ 7 -1CU 2 which determines the be¬ 
havior of the chemistry. The energy balance of high-density 
photodissociation regions is discussed by Burton et al. (1990). 
The main cooling agents of gas at a few thousand Kelvin are H 2 
ro-vibrational lines, followed by H 2 O lines. The heating is pri¬ 
marily provided by H 2 vibrational heating. Far UV photons will 
pump H 2 molecules to a bound excited electronic state (Lyman 
and Werner bands), from which they will fluoresce down to a 
high vibrational level of the ground electronic state. Most of the 
time (85%-90%) the excited H 2 molecule will fluoresce back to 
an low vibrational level while 10-15% of the time the molecule 
will dissociate. At the high densities considered in this study, 
the collision of excited H 2 with other species contribute to the 


heating of the gas and thermalization of the ro-vibrational lev- 


3.2. Chemical model 

We select molecular species composed of five atoms or less 
and containing the elements H, C, O, N, Si, Mg, Fe. The chem¬ 
ical reactions are drawn from the compilation of UMIST rate99 
(Le Teuff et al. 2000). Contrary to the interstellar chemistry in 
cold environments, we expect that reactions with activation bar¬ 
rier, endothermic and three-body reactions dominate the chem¬ 
istry in inner hot disks. Therefore, besides the ion-molecules 
reactions, two new types of reactions are taken into account: 
termolecular reactions, which are catalyzed bimolecular reac¬ 
tions (i.e., 3-body reactions) and collider reactions, which are 
collision-induced dissociations. To emphasize the likely impor¬ 
tance of 3-body reactions, we provide below the dissociation 
sequence of water by a third body which is believed to be the 
initiation reaction in the so-called Bradford sequence. 


H 2 0 + M OH + H (3) 

OH + CO <-> C0 2 + H (4) 

H + H 2 0 <-> OH + H 2 (5) 

H + OH + M —> H 2 0 + M (6) 


where M designates a third body (H, He, H 2 , ...). All together, 
the chemical network consists of 3299 reactions involving 301 
species. Most species have negligible abundance (<10~ 30 w.r.t 
H). No deuterated species are included in our study. 

Another major difference with previous studies of chem¬ 
istry in disks around low- and intermediate-mass young stars 
(Aikawaet al. 1999; Markwicket al. 2002; Nguyen et al. 2002), 
which focus on the chemistry from a few AU (R >10 AU), is 
that the region modeled here is dust-free. The first consequence 
is the absence of the H 2 formation reaction on grain surfaces. 
The second is that far ultraviolet photons are not absorbed by 
dust so that photodissociation is very effective. 

Table 6. Initial abundances. The total oxygen budget is allowed to 
vary to explore the effect of different C/O elemental ratios. Water is 
assumed to be absent at f=0. 


species 

n(species)/n(H) 

H 

1.0 

He 

7.1 x 10' 2 

C 

1.0 x 10~ 4 

N 

2.0 x 10' 5 

O 

1.7 x 10~ 4 

C + 

4.0 x 10' 6 

CO 

8.5 x 10' 5 

h 2 o 

1.0 x IO - 30 
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l uv =1 E8 



l uv =1E7 



l uv =1 E9 



Fig. 5. Relative abundance of various species as a function of time and enhanced interstellar UV field 7uv- The gas is free of dust grains, at 
temperature of 1600 K and has density of 10 12 cnr 3 as deduced from the fitting to the CO and FFO emission lines. CO self-shielding was not 
included in these models. 


3.3. Numerical model 

The model is static and isothermal. The model solves numer¬ 
ically a system of differential equations that simulate the for¬ 
mation and destruction of atomic and molecular species. A 
few solvers of systems of stiff first order differential equations 
were used and their relative speed was checked. The solvers are 
GEAR, VODE (Brown et al. 1989) and two alternative algo¬ 
rithms described by Press et al. (1992). All solvers give similar 
answers for a few test models, ensuring that the results pre¬ 
sented here are not numerical noise caused by a specific solve. 
The adopted solver for the results presented in this paper is 
VODE, the latest incarnation of the backward-differentiation 
method popularized by GEAR. This solver gives the fastest re¬ 
sults on most of the models. We cover a parameters space in 
temperature (1500-5000 K), density (/; h =10' °— 10 14 cm -3 ), UV 
field and CO self-shielding factor. The initial abundance for the 
models is given in table 6. The cosmic-ray flux is set to the in¬ 
terstellar value. Models were run up to one year which is much 
than sufficient for the steady-state to be reached, (see next sec¬ 
tion). 


4. Results of the chemical model calculations and 
discussion 

After presenting the results of the chemical model calculations, 
we will describe the formation route of the main molecular 
species (FE, CO and water) in a hot gas immersed in a strong 
UV field. The possibility that the molecules are formed in the 
outer part of the disk and transported to the inner few AU is 
discussed. Then we consider the effect of CO self-shielding 
against UV photodissociation and non-standard C/O abun¬ 
dance ratio on the molecular abundances. Finally the model 
predictions are compared to the observed abundances. 

4.1. Results 

The abundance of H, FE, O, OH, H 2 O, C, C + , CO, and N are 
shown in Fig. 5 for a gas temperature of 1600 K, gas density 
of 10 12 cm -3 and four values of I uv (10 6 , 10 7 , 10 8 , 10 9 ) and 
no CO nor H 2 self-shielding. First of all, it is clear that steady- 
state is reached in less than a year. This rapidity is reassuring 
because the matter located in the inner few AU of disks rapidly 
falls onto the central star. As a comparison, the rotational pe- 
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riod of the gas at 1 AU around a 6 M 0 star is - 150 days. In 
most cases, steady-state is attained in a few days. The dominant 
species are the atomic H, O, C, C + , N. Molecular hydrogen H 2 , 
CO, OH and H 2 O are the most abundant radicals/molecules. 
Important abundance ratios such as n(H)/n(H 2 ), n(0H)/n(H20) 
and n(CO)/n(HiO) are displayed in Fig. 7 for different com¬ 
binations of h:\- jn\\ and r gas (1600, 2000, 3000 and 4000 K). 
At /uv/hh < 10 -5 , the chemistry is photo-regulated similar 
to a photodissociation region and all newly formed molecules 
are instantaneously destroyed. In the thermal region duv/n H > 
10 5 ), the abundance of CO, OH and H 2 O reduces with increas¬ 
ing gas temperature because, once formed, those molecules re¬ 
act further to synthesize other species. Therefore, the optimal 
temperature for molecular abundance is around 2000+500 K. 


dominates. This reaction has a flat rate of 1.30 x 10 -9 cm 3 s -1 . 
At high temperature, H 2 formation is initiated by the radiative 
association: 

H + + H —> H, + photon (8) 

which has an activation energy of 8500 K, followed by 

Hi + H —> H 2 + H + (9) 

H 2 molecules further react with atomic oxygen or carbon to 
form OH or CH. At temperature higher than 3000 K, H 2 is 
dissociated by the collisions with hydrogen atoms. This latter 
reaction has an activation energy of 55,000 K. At high densities 
(«h >10 12 cm -3 ), H 2 can also be formed by the three-body 
reaction: 


4.2. Hot dense gas chemical network 

The chemical code monitors the main formation and destruc¬ 
tion reactions at designated times. The formation routes for the 
main molecules are discussed below and a schematic is shown 
in Fig. 6. The chemical routes are different for dense and ex¬ 
tremely dense gases. Below densities of around 10 9 cm -3 and 
at high temperature, the chemistry is dominated by two-body 
reactions with activation energy of the order of few thousand 
Kelvin, whereas at very high densities, three-body reactions 
dominate. 



Fig. 6. Schematic representation of the proposed chemical network. 
OH radicals play a central role in this network, v means photodissoci¬ 
ation by ultraviolet photons and M a third body (M=H, He). The main 
routes leading to CO and HiO are shown here. 


4.2.1. Formation of H 2 

Molecular hydrogen is the starting molecule for the chemical 
network. In the absence of dust grains, H 2 is slowly synthe¬ 
sized in the gas phase via three possible route depending on the 
density and temperature. At moderate densities, the reaction 


H + H + H->H 2 + H (10) 

which is a highly exothermic reaction, but the third hydrogen 
atom carries away the excess energy. H 2 reacts swiftly with 
atomic oxygen to from OH which is the key radical for the 
synthesis of CO and H 2 . 

4.2.2. CO formation and destruction 

CO molecules are readily formed for all values of the UV field 
enhancement when the density reaches a value of 10 6 cm" 3 and 
temperatures below 2500 K. The formation of CO is initiated 
by the fast reaction at high temperature (activation energy of 
£ 3=3150 K) 

H 2 +0->OH + H (11) 

followed by the reaction: 

C + OH —> CO + H (12) 

which has a rate of l.lOx 10 -1 °(T/100) 1/2 s -1 . The other chan¬ 
nel 

C + OH —> O + CH (13) 

has an activation energy of 14800 K. 

Above 3000 K, destruction of CO is enhanced because re¬ 
actions with high barrier are activated. Especially, CO is chem¬ 
ically dissociated at gas temperature greater than 3000 K by the 
endothermic reaction: 

H + CO —> C + OH (14) 

This reaction is particularly important since hydrogen is mostly 
in the atomic form. 

At lower UV intensities, CO react with the OH + to form 
HCO + . Here OH + come from the reaction between H j and 
atomic oxygen. 

At extreme densities, CO is formed directly from the 
neutral-neutral addition C+O and the excess energy is carried 
away by a third body, most likely a hydrogen atom. The main 
destruction reaction beside photodissociation is the hydrogena¬ 
tion of CO leading to HCO. If the abundance of H 2 is not suf¬ 
ficiently high 


H - + H -> H 2 + e- 


(7) 


CO + H + M -» HCO + M 


(15) 
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T=1600 K T=2000 K 






Fig. 7. Relative abundance of various species as a function of /uv/uh for gas at temperature 1600 K, 2000 K, 3000 K, 4000 K. 


other possible reactions are 

CO + H 2 + M —> HCO + H (16) 

where M is a third body (likely H or He). Alternatively, with 
low H 2 abundance, CO preferably react with OH to form CO 2 : 

CO + OH —> C0 2 + H (17) 

finally, the reaction: 

HCO + H + M-^H 2 CO + M (18) 

convert HCO into H 2 CO. Therefore, HCO and H 2 CO have non- 
negligible abundance. 

4.2.3. Water molecules 

The chemical reactions leading to water molecules in cold and 
hot regions are different. In quiescent molecular cloud, the gas 
temperature is below 230 K and gas-phase chemistry gives typ¬ 
ical H 2 O abundances of a few x 10~ 7 (e.g., Lee et al. 1996; 
Le Teuff et al. 2000). Water is formed by a sequence of ion- 
molecules reactions starting with O + H | or 0 + + H 2 lead¬ 
ing to OH + . Then rapid H-abstraction reactions with H 2 pro¬ 
duce H 3 0 + which dissociatively recombines (with electrons) 


to H 2 O. Another potential source of gas phase water molecules 
is the evaporation of icy grains in regions where T >100 K. 
The reaction scheme in cold regions has been proven to be 
too effective in producing H 2 O when compared to observa¬ 
tions of high-mass star forming regions by Snell et al. (2000). 
At high temperature (7„ as >230 K) and densities lower than 
10 15 cm -3 , water molecules are predominately formed via the 
radical-molecule reaction: 

H 2 +0H-^H 2 0 + H (19) 

since OH is very abundant. The importance of this reaction 
has already been realized for water formation in shocks (e.g., 
Elitzur 1979; Bergin et al. 1998). Water molecules are mostly 
destroyed by photodissociation and by reaction with C + : 

H 2 0 + C + -> HOC + /HCO + + H (20) 

The Bradford sequence for the synthesis of water still plays 
a small role compared to the radical-molecule formation route 
described above at the densities attained in inner disks. The 
reaction network (CO and H 2 O) discussed here and summa¬ 
rized in Fig. 6 is similar to those found for post-/-shock chem¬ 
istry (Hollenbach & McKee 1989). It may therefore be diffi¬ 
cult to distinguish between a post-shock or a quiescent high- 
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temperature chemistry origin for many molecules. The produc¬ 
tion of water vapor by sublimation of icy mantle when the dust 
temperature exceeds 100 K is often advocated to explain large 
quantity of FFO in the gas phase. 



Fig.9. 1 FO/CO versus FFO/OH diagram for various 7 uv/wh and gas 
temperature. 


4.3. Endogenous versus exogenous formation of 
molecules 

Two main scenarii can explain the presence of molecules in the 
inner disks around young massive stars. The first model pro¬ 
poses that the molecules are formed at large distances from the 
central object and subsequently brought to its vicinity by ac¬ 
cretion. We call this possibility the exogenous scenario. In the 
exogenous model, molecular species can also originate from 
the evaporation of icy grain mantle. The second model states 
that molecules can be readily synthesized locally in the inner 
disks. This is the endogenous model. All models reach steady- 
state in at most one year and the effects of the initial molecular 
abundances are minimal. 

4.4. Effect of CO self-shielding 

At large column densities, CO molecules can self-shield from 
photodissociation (e.g., van Dishoeck & Black 1988). The 
nominal column of CO is 10 15 cm -2 for a turbulent width of 
5 km s 1 but this value was derived for gas at low temperature 
and may not be valid for gas above 1500 K, because highly 
vibrational excited absorption levels were not taken into ac¬ 
count in models of CO self-shielding. The nominal column 
of CO for self-shielding varies as 1/Av. For a very turbulent 
gas with Av= 30 km s~ 1 , the nominal column of CO is in¬ 
creased by a factor ~6. Nevertheless, we ran models with in¬ 
creasing CO column density N( CO) to account for probable in¬ 
crease of the value of the nominal CO for self-shielding needed 
to simulate the condition in inner disks. The other parame¬ 
ters are 'f gas = 1600 K, » H =10 12 cm -3 and 7uv=10 8 . Figure 10 
shows that above some threshold value (A(CO)>10 15 cm -2 ). 



Fig. 10. Effect of the CO self-shielding on the relative abundances. 
The self-shielding factor depends on the CO column density N{ CO). 
The parameters of the model are r gas =1600 K, n H =10 12 cnr 3 and 

/uv=10 8 . 



Fig. 11. Effect of the oxygen elemental abundance on the relative 
abundance of CO, OH and H20. X D (0) is the nominal value of 
the Oxygen elemental abundance. The parameters of the model are 
r gas =1600 K, n H =10 12 cm- 3 and / uv =10 s . 

CO molecules are self-protected against photodissociation and 
its abundance can reach the cold quiescent molecular cloud 
value of ~10 -4 . Surprisingly, OH and H 2 O abundances seem 
to decrease when CO abundance rises steeply. The H 2 O/CO ra¬ 
tio drops dramatically from few 10 4 to few 10 7 . From Fig. 8, 
it is clear that the effect of CO self-shielding is most notice¬ 
able for /uv/»h > 10 3 (the photodominated regime). The CO 
abundance reaches similar value than in the thermal region, i.e. 
where 7jjv/»h < 10 -5 . 

4.5. Effect of C/O abundance ratio 

By pushing the dust grains beyond the dust sublimation radius, 
the radiation pressure will create a separation between gas and 
dust. The segregation of the gas and dust results in an increase 
of the elemental C/O ratio because the silicate grains, which 
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Fig. 8. Effect of CO self-shielding on abundance ratios. The left panel shows the addition of a strong self-shielding on the chemistry while 
the right panel is the same than the upper right panel in Fig. 7. As expected the changes happen only on the HiO/CO and CO/H ratios. The 
observed value of FFO/CO (~ 5.4 10~ 4 ) is indicated by a grey strip. 


contain 10-20% of the available oxygen, do not sublimate. The 
effect of an increase in the elemental C/O ratio above the solar 
value of 0.4 would be to lower the water and molecular oxy¬ 
gen abundance. In principle, the elemental abundances can be 
deduced from the study of photospheric lines of the star if the 
star is not too embedded. Terzieva & Herbst (1998) discussed 
the different effects on the chemistry caused by variations in 
the C/O elemental abundance but their study is limited to low 
temperature chemistry where elements can be depleted onto 
grains. Here we allow the oxygen elemental abundance to vary 
from 0.1 to 1.0 times the nominal value listed in Table 6. The 
model results at steady-state are shown in Fig. 11. As expected, 
the relative abundance of CO, OH and H 2 0 diminish with de¬ 
creasing oxygen elemental abundance. The rate of decrease is 
similar for the three species. Therefore, the relative abundances 
HtO/CO and FFO/OH are not significantly affected by chang¬ 
ing the oxygen abundance. The simultaneous decrease in CO 
and HiO abundance with lower oxygen elemental abundance 
can be ascribed to the formation route of both molecules with 
are formed via the radial OH. This is in contrast with the chem¬ 


istry that occurs in quiescent molecular clouds where H 2 0 is 
synthesized at a much lower rate than CO, which results in a 
higher CO/H 2 0 ratio as the oxygen elemental abundance de¬ 
creases. 

4.6. Comparison with the observations 

A large range of gas temperature (1600-4800 K) has been in¬ 
ferred from the modeling of CO bandhead emission toward 
four high-mass young stars (Bik & Thi 2004). The abundance 
of CO starts to decrease at temperature above 2500 K owing to 
collisional dissociation of CO and to the reaction of CO with 
OH to form C0 2 . 

The observed H 2 0/CO ratio in 08576nr292 is 5.4 x 10" 4 . If 
CO self-shielding is taken into account, the value for /uv/«h is 
not well determined because the theoretical curve for H 2 0/CO 
is relatively flat around few times 10 4 (upper left panel in 
Fig. 7). The CO abundance is 10 7 —10^ 4 and H 2 0 abundance 
is 10~ n -10 -8 . Interestingly, when CO is allowed to self-shield, 
we can better infer 7 uv/«h ~10 7 for a gas at T=1600 K from 
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the observed value of FFO/CO since the "degeneracy” in the 
H 2 O/CO curve in the right panel of Fig. 8 is lifted when CO 
self-shielding is taken into account in the modeling. 

The absolute abundances are low in comparison to that 
in the envelope of high-mass young stellar objects (Boonman 
et al. 2003). From the CO abundance in the unshielded case, we 
infer a total hydrogen column density of 3.9 x 10 25 —10 28 cm 2 
or surface density So = 6.5-65,000 g cm -2 . In the case CO self¬ 
shielding effects are included in the modeling, CO/H reaches 
the maximum relative abundance of 10 4 and the total hydro¬ 
gen column density amounts to 3.9 x 10 25 cm 2 or a surface 
density of So = 6.5 g cm 2 . We provide a large possible range 
for the derived total hydrogen column density because current 
CO self-shielding models may not be applicable to the situation 
encountered in inner disks. However, it is noteworthy that the 
upper range of column density (V(H)>10 26 cm 2 ) would imply 
large opacity in the near-infrared due to H , which would pre¬ 
vent the detection of optically thin H 2 O lines. The surface den¬ 
sity is comparable to the values adopted in disk models around 
T Tauri and Herbig Ae stars which is So —10 2 -10 4 g cm 2 at 
1 AU. 

We summarize the chemical modeling in Fig. 9. This fig¬ 
ure shows the H 2 O/CO versus H 2 O/OH abundance ratio. The 
figure provides a simple way to deduce physical parameters 
such as Tgas and /uv/«h from CO, OH and H 2 O abundance. 
Unfortunately, OH lines were not covered in the K -band spec¬ 
trum of 08576nr292. 

4. 7. Comparison with the inner disk chemistry around 
low-mass young stars 

Glassgold et al. (2004) studied the thermal-chemical structure 
of the inner disks around young low mass stars. Although in 
low-mass young stars X-rays drive the chemistry rather than 
the UV, they found that the disk atmosphere may contain a sig¬ 
nificant amount of warm ( T =500-2000 K) CO gas similar to 
our results. This similarity is not surprising since at the high 
temperatures and densities encountered in inner disks around 
low- and high-mass young stars, the chemistry is dominated by 
the same neutral-neutral reactions. At very high densities, the 
UV and X-ray photons provide the energy needed to maintain a 
high temperature. The photodissociation by X-ray or UV pho¬ 
tons of the molecules are compensated by rapid formation in 
dust-free regions. In a dust-free region X-ray and UV photons 
are absorbed by gas while in a dust-rich region UV radiation is 
rapidely attenuated by dust absorption while the X-ray can pen¬ 
etrate deeper to influence the ionization fraction (e.g., Semenov 
et al. 2003). 

5. Conclusion and future prospects 

We have shown that hot-water vapor emission is present in the 
close vicinity (R <5 AU) of the high-mass young stellar object 
08576nr292. The water vapor as well as the hot CO are likely 
located in a disk in Keplerian rotation. 

The detection of optically thin CO and H 2 O lines at 1500- 
1600 K suggests that dust grains are absent in that tempera¬ 
ture range in the disk. This may stem from the fact that dust 


grains are pushed further out by radiation pressure or that most 
dust has evaporated. The region probed by hot CO and H 2 O is 
limited to the inner few AU of the disk. Alternatively, the gas 
and dust may be thermaly decoupled with the gas being much 
warmer than the dust. 

To probe cooler gas, which is located between 1 and ~50 
AU, fundamental transitions of CO and lower J lines of H 2 O 
appear to be more appropriate although the line over continuum 
contrast would be smaller owing to the warm dust emission. 
Nevertheless, hot lines allow the gas-rich and probably dust- 
poor region to be investigated. Higher resolution spectrometers 
(R- 100,000 or Av=3 km s _1 ) such as CRIRES planned for the 
VLT. will resolve the individual H 2 O line profiles. The combi¬ 
nation of medium resolution (R- 10,000) and high spatial res¬ 
olution achieved by AMBER will help to determine the exact 
location of both the CO and H 2 O emissions. 

Chemical models have been developed for the dense and 
hot inner disk around young high-mass stars. The chemistry 
is dominated by neutral-neutral reactions and steady-state is 
rapidly attained. H 2 formation is hampered by the lack of dust 
grains and reach relatively low abundance. The radical OH is 
central to the chemical network playing the same role than 
HI for cold chemistry. The molecular abundances depend in 
general on the gas temperature and the UV flux over den¬ 
sity ratio and show a large range of values. Photodissociation 
affects more H 2 , OH and H 2 O than CO because the latter 
species can self-shield very efficiently. The observed H 2 O/CO 
can be reproduced at the temperature of the gas (1600 K) and 
7 uv/«h -10 Observations of OH are warranted to allow bet¬ 
ter constraints on the UV radiation field. Future work on the 
modeling side includes the combination of the chemical model 
presented here with a dynamic disk model around young high- 
mass stars. 

The detection of molecules such as CO and H 2 O at high 
temperature in the close vicinity of a young massive young 
stellar objects opens the possibility to study the physics and 
chemistry in the region where matter from the disk falls onto 
the star. 
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